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Abstract 

The effects of a non-magnetic Zn impurity substituting an in-plane Cu are studied by solving the 
Bogoliubov-de Gennes equation self-consistently which is derived from the t-t'-J-U Hamiltonian 
with all the allowed order parameters included. The Zn impurity, modeled in terms of a potential 
scatterer in unitary limit, induces local staggered magnetic moments around itself, and the cal- 
culated NMR shifts from the induced moments are in agreement with the experimental Cu NMR 
spectra. We also note that the experimentally observed negative slope of the tunneling conductance 
can result from the next-nearest hopping t' . 
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After more than a decade of intensive research on the high T c superconductors, we are still 
far from a coherent understanding of the underlying physics |L|. One informative route of 
study is provided by the impurities like Zn and Ni, which substitute for Cu sites in the Cu0 2 
planes where the superconductivity (SC) is believed to emerge from. Therefore, the impurity 
induced information about the SC and magnetism and their correlation can provide a clue 
to the underlying physics of the high T c SC in that a common feature of the cuprates is the 
proximity of the antiferromagnetism and SC. The high resolution atomic scale probes on the 
impurity substituted cuprates provide important information in this regard. For instance, 
Pan et al. reported the scanning tunneling microscopy (STM) study on a Zn substituted 
Bi 2 Sr 2 CaCu 2 8 +<5, that is, the tunneling conductance as a function of the atomic position 
and bias voltage, G(f, V) ||. They observed among others that (a) at r = Zn site, there is 
a strong peak in G(f, V) as a function of V near zero bias voltage (V ~ —1.5 meV), (b) at 
V = —1.5 meV, the real space distribution of the intensity G(f, V) shows local maxima at 
next-nearest Cu sites from the Zn, and (c) the spatially averaged conductance / drG(r, V) 
has a negative slope as a function of V. Also, several groups reported NMR experiments on 
Zn || |j, Li ||, or Ni || substituted YB 2 Cus07_ :r . It was found that (d) local staggered 
moments were induced and SC was destroyed around the impurity site, and away from there 
the induced moments become vanishing and the pairing amplitudes heal back to the bulk 
value. 

Interpretation of these observations at present is controversial. For the observations 
(a) and (b), some groups proposed the blocking or interference effects |?J |§, but other 
groups introduced the Kondo effect due to the impurity induced magnetic moments ||. 
For the observation (c), Rantner and Wen suggested from a £77(2) slave boson mean- field 
calculation that the asymmetric tunneling background may be due to the strong onsite 



repulsion fllOfl . Although there are quite a few theoretical works on the magnetic or non- 
magnetic impurity effects on superconductivity |Tl|, [12|, 13|, |14 |, no theoretical calculation 
of the induced magnetic moments due to a non-magnetic impurity in a SC state has been 
reported yet. We report such a calculation in this paper. We study the t-t'-J-U Hamiltonian 
with the Zn impurity modeled as a potential scatterer near the unitary limit to understand 
the observations (c) and (d) listed above. 

In order to address the observations (c) and (d), we consider the t-t'-J-U model given as 
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follows: 

<i,j>a ' <i,j> V 

+ UY,ni\nil + VtmpY, c t a c io<T, (!) 

i cr 

where n ia = cf a c ia , Hi = J2a^ia, Si the spin- 1/2 operator, /i the chemical potential, U the 
onsite Hubbard repulsion, Vi mp the impurity potential at the Zn impurity site z'q, and J is 
the nearest neighbor exchange integral. The <> and <c;§> stand for the summations up to, 
respectively, the nearest and next-nearest neighbors, and is equal to t for the nearest and 
to t' for the next-nearest neighbors. The normal state band structure is then given by 

6^ = — 2t[cos(fc x a) + cos(k y a)} — At' cos(k x a) cos(k y a). (2) 

Note that we included both U and J in the Hamiltonian. This was motivated by the recent 
work by Daul et al. who suggested that the t-U-J model be more appropriate than the t-J 
or Hubbard models to describe the real cuprate system by allowing the charge fluctuations 
with the finite Hubbard U |15| . 



The t-t'-J-U model of Eq. ([!]) was treated with the Hartree-Fock-Bogoliubov (HFB) ap- 
proximation to yield the Bogoliubov-de Gennes (BdG) equation in the real space. It was 
then solved via iterations with all the allowed order parameters determined self-consistently, 
that is, the (i-wave pairing, bond, charge, and spin order parameters. We found that (1) local 
antiferromagnetic moments are induced due to a non-magnetic impurity, (2) the tunneling 
conductance shows a negative slope as a function of bias voltage V due to the next-nearest 
neighbor hopping t', and (3) the local density of states (LDOS) at the 4 nearest neighbors 
shows a peak near zero energy. The calculated NMR spectra from the distribution of the 
induced magnetic moments is found to be in agreement with the experimentally observed 
63 Cu NMR spectra M. There is some controversy in understanding the experimental obser- 
vations (a) and (b) with the result (3) as mentioned before. Here, we will report the results 
(1) and (2) in more detail to understand the observations (c) and (d). 

The BdG equation derived from the Hamiltonian of Eq. (|l|) within the HFB approxima- 
tion can be written as 
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where 



and the subscripts r and r' are, respectively, the nearest and next-nearest neighbors. The 
order parameters, <i-wave pairing A^-, bond-order Xijai s Pi n U P an d down charge densities 
n ia , were determined from the following self-consistency conditions: 

A tj = J({Ci^Cji) + (CjiCil)) = jY^[Uj, n Vi,nf(-E n ) +U i>n Vj )n f(-E n )], 

n 

Xitf = (4^n) = J2 U i,n U j,nf( E n), 
n 

Xijl = (ctl C jl) = Y, V tnVj,nf(-E n ), 
n 

= (n iT ) = J2 U tn U i,nf(En), 

n 

n H = (nil) = J2 V tn V i,nf(-En), (5) 
n 

where (• • ■) represents the thermodynamic average and f(oS) = 1/[1 + exp(/3u>)] is the Fermi 
distribution function. Charge and spin densities are then given by rii = + and 
i(ra iT - Tin), respectively. Eq. (H) with the self-consistency conditions of Eq. (^) was 
solved via iterations with a unit cell of 20 x 20, or 30 x 30 sites using a periodic boundary 
condition with 5x5 wave-vector points in the first quadrant of the Brilouin zone. The 
LDOS which is proportional to the tunneling conductance, and local magnetic moments of 
the z-th site can be obtained from 

D(i,u)=^2 [\u^ nk \ 2 5(uj - E nk ) + \v i>nk \ 2 6(uj + E nk ) , 

n,k 

S zfy=2 ^ [\ U iM\ 2 f( E nk) ~ \Vi,nk\ 2 f(-E nk ) ■ (6) 
n, k 

A few remarks about the Hamiltonian are in order. The bond-order, although not essen- 
tial, was included because (a) it is an allowed order parameter in the t-t'-J-U model unlike 
in a model with the onsite repulsion U and next-nearest-neighbor attraction V (t-t'-U-V , 
or minimal model), and (b) it was found important for stablizing the bond-centered stripes 
where the more sophisticated methods like the density matrix renormalization group predict 



so and the half-filled metallic stripes as observed experimentally [17]. We wished to use a 



model of more general validity applicable for both the impurity and stripe problems. This 
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FIG. 1: The tunneling conductances in normal (dashed) and SC (solid curve) states. In the 
inset, TDOS in SC state in the main figure is shown in more detail around the gap feature and 
is compared with the Pan et al. experiment (dashed curve) and Rantner and Wen calculation 
(dot-dashed). 

is the main reason we employed the t-t'-J-U model in the present study. On the other hand, 
the physics discussed here may also be described in terms of the minimal model because the 
bond-order is not essential for the moment induction. 

Let us first consider the asymmetric tunneling background of the observation (c). We 
take the parameters t = 0.4 eV, t' = —0.45 t, U — 3 t, and J = 0.75 t, which are appropriate 
for E^S^CaC^Os+s- Rather a large value of J was necessary to fit the pairing amplitude 
A = 48 meV observed in the Pan et al. STM experiment. To illustrate the point in a 
simple way, we put S z = 0. Fig. 1 shows the tunneling density of states (TDOS) in normal 
(dashed) and SC state (solid line) at the 16 % doping concentration. The two-dimensional 
tight-binding band given by Eq. (g) has a van Hove singularity (vHS) due to a logarithmic 
divergence. The vHS, or a peak in DOS occurs at the band center for t' = 0. But, for t' < 
the peak in the DOS moves to a lower energy relative to the Fermi level and, consequently, 
the DOS has a negative slope around the Fermi level {u = 0) as a function of the bias 
voltage as shown in the Fig. 1. Small offset of the vHS peaks between the normal and 
SC states are due to the pairing and the difference in the bond orders, and the V-shape 
around u = of the solid curve is due to the <i-wave pairing. This is shown in more detail 
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around the gap feature in the inset. The solid, dashed, and dot-dashed curves represent, 
respectively, the present calculation, Pan et al. data, and Rantner and Wen calculation. The 
agreement of the present calculation with the Pan et al. data is quite satisfactory without 
fine adjustments. It seems natural, therefore, to understand the observed negative slope of 
the tunneling background O, [TB[ in terms of the particle-hole symmetry breaking around 



the band center because of the next-nearest hopping integral t'. Rantner and Wen, on the 
other hand, argued that the asymmetric tunneling background can be explained in terms 



of SU(2) slave-boson approach [10]. But the present explanation is more straightforward, 
and fits the experimental observation much better than their calculation. Interactions may 
bring in some corrections in the tunneling conductance but the background slope will be 
mainly set by the particle- hole symmetry breaking due to the t' . We therefore believe that 
the asymmetric tunneling background is most probably due to the next-nearest hopping 
integral t' . 

Now we turn to the main result of the present work, that is, the magnetic moment forma- 
tion due to a Zn impurity substituting an in-plane Cu site of the cuprate superconductors. 
We took J = t, T = 0.0005 t, and all other parameter values the same as in Fig. 1. A 
slightly larger value of J was taken to show the doping induced features more clearly. This 
set of parameters yields the pairing A = 0.312 t, magnetic moment S z (i) = 0, and bond 
order Xa= 0.18 t. The doping level was set such that (i) the self-consistently determined S z 
vanishes without the impurity and (ii) the chemical potential is higher than the vHS energy. 
The requirement of (i) ensures that the moments do not exist without the Zn impurity, 
and (ii) ensures that the tunneling conductance background has a negative slope around 
the Fermi level. In the present HFB mean-field calculation, these two requirements are met 
around 40 — 55 % doping for t' = —0.45 t depending on U and J. This extremely high dop- 
ing concentration, an artifact of the mean-field calculation which could be cured by going 
beyond the mean-field theory, was caused by the requirements of negative tunneling con- 
ductance slope and the large (realistic) t' . These constraints were not properly considered 
previously. Although the present calculations can also be done around 15 % doping without 
the constraints, they were done at a high doping level with the experimental constraints 
satisfied. The reason is that the induced moments are not affected much by the largeness of 
the doping level since they are mainly determined by the competition of the electron kinetic 
energy and the effective exchange integral. This has been checked by direct numerical calcu- 



6 



30 



FIG. 2: The staggered magnetic moments S z (i) induced around an non-magnetic impurity at 
T = 0.0005 t, for J = t, t' = —0.45 t, and U = 3 t. Note that the induced moments reside mainly 
on the 4 nearest neighbors and exhibit a Friedel-like oscillation. 

lations without the requirements, and is further strengthened by the satisfactory agreement 
of the computed NMR spectra from the induced moments with experiments as is discussed 
in Fig. 3. On the other hand, the resonant energy and peak strength depend on the param- 
eters of the model like the t', Vi mp , and doping concentration somewhat sensitively as noted 
previously |]I5fl , and some care must be taken before comparing different calculations. 



We took Vimp = —100 t to model a Zn impurity placed at i — (15,15) in a 30 x 30 unit 
cell. Fig. 2 shows the staggered magnetic moments S z (j) = (— l)^ x+lyS) S z (i) induced by the 
Zn. Antiferromagnetic moments develop with the peaks at the four nearest sites around the 
impurity site. The total spin J2i S Z (J) = \h. The magnetic moment vanishes at the impurity 
site, and the Zn remains non-magnetic. The Friedel-like concentric oscillation with a period 
of 3 — 4 lattice constants can clearly be seen. The oscillation exhibits an angular anisotropy 
as expected; along the anti-nodal direction it gets more suppressed because of the pairing 
gap and along the nodal direction it shows more pronounced oscillation. The d-wave pairing 
gap, A d (i) = [Afj +X + A^_ x - A^ +y - A^_ y ]/2, remains 0.312 t for most of unit cell, but 
on the impurity site A^(ij vanishes. The superconductivity is severely destroyed by the 
non-magnetic impurity within about 2 lattice constants around the impurity site. 

From the Zn impurity-induced S z (i) shown in Fig. 2, we have also calculated the broad- 
ened NMR spectra which may be compared with the experimental 63 Cu NMR spectra. 
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FIG. 3: The Cu NMR spectra calculated from the distribution of the induced magnetic moments 
due to a Zn impurity as shown in Fig. 2. 

The NMR spectra is just a histogram of the frequency shifts given by hi = A°™ stte S ' z (i) + 
BJ2 T S z {i + f), where A° a nsite = -13.8 T/fi B and B = 3.45 T/fi B for the magnetic field 
H || c-axis The T-dependence of the induced moments were calculated for several dif- 
ferent T with a 20 x 20 unit cell. The distribution of the induced moments in 20 x 20 
unit cell is almost the same with the 30 x 30 one because the magnetic correlation length 
is shorter than the unit cell size. The results at T = 0.1 T c , 0.2 T c , 0.5 T c , and 1.5 T c are 
shown in Fig. 3. The width of the NMR spectra gets narrower as T is increased because the 
induced moments are decreased. Compared with the NMR spectra [Fig. 1(a) in Ref. ||]], 
the results are quite satisfactory, which gives a reliability to the present calculation of the 
induced magnetic moments. 

There are two factors that contribute to the induction of the magnetic moments by a 
non-magnetic impurity: the localization of an electron and destruction of superconductivity 
around the impurity site due to the deep impurity potential of the Zn. The magnetic moment 
in a normal state is determined by the competition between the effective exchange integral of 
2J + U and the electron hopping motion of t. When an electron is localized by the impurity 
potential, the electron hopping gets suppressed, which favors the appearance of magnetic 
moments around the impurity site, as was studied by Odashima and Matsumoto in a normal 



state [19|]. In a SC state, the appearance of magnetic moments can also be suppressed by 



the onset of superconductivity. Therefore, when the superconductivity is destroyed by an 
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impurity, magnetic moments can appear around the impurity site. Both of above two factors 
contribute to the appearance of the magnetic moments induced by the Zn impurity in a SC 
state. Similar ideas were suggested by various groups in the context of spin gaped systems 
PH [?T| ]. The idea was that the freed spins created by an impurity abide by the impurity 
site because of the confinement provided by antiferromagnetic correlations and show up 
around the impurity. The present theory can readily be applied to the magnetic impurity 



effects such as the recent STM experiment on Ni substituted Bi 2 Sr 2 CaCu208+5p2|. We are 
currently carrying out such calculations and will report the results elsewhere. 

To summarize, we considered the t-t'-J-U Hamiltonian to study the non-magnetic Zn 
impurity effects in d-wave superconductivity. The Bogoliubov-de Gennes equation in the 
real space was derived from the Hamiltonian using the Hartree-Fock-Bogoliubov approx- 
imation, and was solved self-consitently with the d-wave pairing, bond, charge, and spin 
order parameters included. It was demonstrated that a non-magnetic impurity induces the 
antiferromagnetic moments around itself in a <i-wave superconductive state. The calculated 
NMR spectra due to the induced moments are found to be in good agreement with the 
experiments. It was also proposed that the negative slope of the tunneling conductance may 
most probably come from the next-nearest hopping integral t' . 
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